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We analyze electron paths in a solid-state double-slit interferometer based on the 
two-dimensional electron gas and their mapping by the scanning gate microscopy 
(SGM). A device with a quantum point source contact of a split exit and a drain 
contact used for electron detection is considered. We study the SGM maps of source- 
drain conductance (G) as functions of the probe position and find that for a narrow 
drain the classical electron paths are clearly resolved but without any trace of the 
double-slit interference. The latter is present in the SGM maps of backscattering 
(R) probability only. The double-slit interference is found in the G maps for a wider 
drain contact but at the expense of the loss of information on the electron trajectories. 

Stability of G and R maps versus the geometry parameters of the scattering device 
is also discussed. We discuss the interplay of the Young interference and interference 
effects between various electron paths introduced by the tip and the electron detector. 

I. INTRODUCTION 

The scanning gate microscopy (SGM) is an experimental technique 1 ^ that probes transport 
properties of devices with the two-dimensional electron gas - buried shallow beneath the surface 
of the sample - by the charged tip of the atomic force microscope. The technique has been 
used in particular for investigation of quantum point contacts - the most elementary quantum 
transport device^®. The conductance maps gathered with the SGM contain characteristic 
oscillation of the period of half the Fermi wavelength that appear due to the interference 
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of electron wave functions incoming from the quantum point contact and backscattered by 
the tip 5 13 . In our recent papei ! 14 we have proposed a system with split source channel for 
observation of the double-slit interference. The Young interference should be present in the 
SGM maps provided that the transport in the channel which feeds the split source occurs 
in the lowest subband of lateral the quantization 14 . The proposal of Ref . 14 and most of the 
experimental studies 5 8 dealt with systems in which the electron after passing through the 
constriction defining the quantum contact enters an infinite half-plane. In this work we consider 
the possibility of observation of electron paths in the context of the double-slit interference. 
According to the which-path thought experiment by Feynmaii 15; determination of the slit the 
electron goes through destroys the double-slit interference. Here, we demonstrate that SGM can 
be used for detection of classical electron paths - although without indicating the one taken 
by the electron - with a simultaneous resolution of the double slit interference effects. The 
SGM was early used for detection of the semi-classical electron trajectories as deflected by the 
Lorentz forcd- 16 , including observation of the skipping orbit ^ 7 18 for systems with an additional 
confinement. Observation of electron paths requires placing an electron detector in the system 
- usually another QPC serving as the drain channel 16 18 . In the present paper we consider a 
system with a split QPC serving as the electron source and the second QPC used for detection 
of the electron passage. We demonstrate that the source-drain conductance maps for a narrow 
drain detector indicate the classical paths but miss the double-slit interference. For enlarged 
drain width the double-slit interference appears in the map but the image of the paths is lost. 
The simultaneous observation of both the semi-classical paths and the Young interference is 
possible when besides the map of source-drain conductance one considers the SGM maps for 
backscattered electrons - or equivalently - the maps for conductance between the source and 
the rest of the system excluding the drain detector. We also discuss the electron paths in terms 
of the quantum consequences in various interference effects. In particular we also point out 
that backscattering by the tip induces interference between the slits with a characteristic SGM 
pattern that is superimposed on the Young interference image. 
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FIG. 1: a) Sketch of the system considered in this paper with the double-slit system - source channel 

- QPC1 of width tCqpd that splits the incoming electron wave function into two beams. The width 
of the QPC1 slits equal to the width of the input channel is kept at tc qP ci = 40 nm all along this 
work. The slits of QPC1 in this Figure are separated a distance of d qpc i = 160 nm. In the rest of 
the paper d qpc i = 400 nm is mainly used. The outgoing current is gathered by the drain channel 

- QPC2 at right of width w qpc 2 to determine source-drain conductance G. Additionally the current 
between the source and the rest of the system excluding the drain detector is measured for evaluation 
GVest leakage conductance. The sample is assumed large and the blue dashed line indicates the ends of 
the computational box where transparent boundary conditions are applied. The distance between the 
computational box and the edge of the sample (Id) is assumed much bigger than the coherence length 
Id ^ A C oherence- The blue rectangle, between the left and right QPCs, shows the area of calculated 
SGM conductance maps. The distance between QPC1 and QPC2 is 600 nm. b) The absolute value of 
the scattering wave function \i/j\. c) Real part of ip. d) Real part of the wave function calculated from 
equation (J9| as superposition of three point sources in the centers of the slits. 


II. MODEL 


We consider an experimental situation that is depicted in Fig. [l](a) . The current is fed by 
the source contact to the channel that is filtered by the split quantum point contact QPC1. The 
second QPC (QPC2) serves as an entrance to the drain contact at the right side of the figure 
that will be used for electron detection as in Refs. 16 18 l or more recently in Ref. 19 . The effective 
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width of the source channel w qpc \ and the width of both slits at the QPC1 side is taken equal 
40 nm, which for the considered Fermi energy transmits the current in the lowest subband only. 
According to the previous papei 14 the lowest subband transport on the input part is necessary 
for the Young interference to be observed. For a larger number of incident subbands the 
conductance map for the double slit system is a simple sum of maps for separate QPCs 1 ^ since 
the double-slit interference disappears in the Landauer summation over the incident subbands. 
The electron leaving QPC1 enters a finite but large sample filled by 2DEG. By large we mean 
that the G distance between the region of interest for the SGM and the edges of the sample is 
larger than the coherence length, so that the interference with the edges can be neglected. The 
source-drain voltage is considered low enough for the linear transport conditions to occur. The 
current passing to the drain is measured in order to evaluate the G conductance. The sample is 
grounded by a large reflectionless contact and this current is measured to evaluate the leakage 
conductance G rest . From conservation of the current and for M in subband in the input channel 
(we consider mostly M m = 1) we have 

Min = R + G + G res t, (1) 

where R = ^ P( IS , where Pb s is the backscattering probability. Thus, R can be determined 
when G and G rest are measured. In the following we discuss numerical results for R and G. 

In the calculations we focus on the region marked by the dashed lines in Fig. [IJa). We 
assume that the transport is coherent within the computational box. At the dashed lines we 
apply transparent boundary conditions, so that the region of interest is effectively open, in 
contrast to systems with a pair of QPCs used as source and drain for a closed stadium that 
are studied by SGM in e.g. Refs. 20 21 . In this paper we set the distance between the input 
slits dqpci = 400nm, which is large enough to distinguish trajectories of the electron arriving 
to QPC2 from one of the input slits. In Fig. [l] a smaller value of d qpcl = 160nm is used for 
illustration. The total computational box covers the region where the transparent boundary 
conditions are applied of size 800nmxl000nm and the region where the scans by the SGM are 
taken is significantly smaller: 400nm x 600nm - see Fig. [l](a). 

In order to simulate the propagation of the Fermi level electrons within the region inside the 
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computational box (dashed lines in Fig. [l](a) ) we consider coherent transport as described by 
the effective-mass Schrodinger equation, 

{~2^ v2 + yti = E F'l’(x,y), (2) 

where m e g = 0.067mo is GaAs electron effective mass, Ep is the Fermi level energy, and potential 


Kip (x,y) 


_Gnp_ 

1 + [(z - ^tip) 2 + (y- 2 /ti P ) 2 ] /dtip 


( 3 ) 


describes the effective perturbation induced by the AFM tip with amplitude f7 t i P = 15 meV and 
width d tip , localized above point x tip , y tip . This type of effective Lorentzian-shaped perturba¬ 
tion, which results from the screening of the charge on the tip by the 2DEG electron gas located 
under the sample surface, was obtained previously in our self-consistent Schrodinger-Poisson 
calculations 22 23 . The width of the potential is of the order of the distance between the tip and 
the 2DEG 23 . In the SGM experiments the 2DEG is buried at least 25 nm below the surfaced, 
and the minimal distance of the tip to the surface applied in SGM is 20 nmP, hence a minimal 
realistic value of du p is about 50 nm. In this paper we consider two widths of the tip: a small 
one d t i P = 10 nm which is useful for the initial discussion since it sets the precision of the 
determination of electron paths, and a realistic one - du p = 50 nm. 

For simplicity in order to define the contacts we assume hard-wall boundary conditions 
on QPC1 and QPC2. In experimental setups the QPCs are usually defined electrostatically 
by potential applied between the gates so that the QPC potential has a saddle point profile. 
Nevertheless, we do not discuss conductance quantization as function of the QPC width and 
for a fixed number of transmitting subbands the hard-wall potential gives qualitatively similar 
results to a saddle point potential which was applied in a previous work 2 '. The transparent 
boundary conditions at the blue dashed line of Fig. [lja)) were introduced by the method 
described in Ref. 2 '. 

We work within the Landauer approach for a zero temperature in the finite difference 
implementation 28 of the quantum transmitting boundary method 29 30 . Within the input chan¬ 
nel - before the splitting ending with the two QPC1 slits, the wave function at the Fermi level 
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is given by the superposition of incoming and backscattered transverse modes 

M qpcl 

VW(^2/) = Y { a ke lkx xT cl (y) + he~ lkx x qP k 1 (y)} 

k= 1 
+00 

+ £ h^’xr'W; ( 4 ) 

k=M qpcl+1 

where the last term corresponds to the summation over the evanescent rnodetP®- anc [ M qpcl 
is the number of current propagating transverse modes y)P cl in QPC1. We consider a single 
subband in the QPC1 channel and an arbitrary number of subbands in the QPC2 channel. 
The summation runs over the Fermi level wave vectors k at subsequent lateral subbands. The 
transport from the input channel to the two-slits is non-adiabatic with a pronounced backscat- 
tering. Note, that the central island [Fig. 1(a)] used as a beam splitter in the presence of a 
strong external magnetic held is likely to form a quantum Hall interferometer as discussed in 
Refs.™ 

For QPC2 we assume the boundary conditions of form 

M qp c2 

^qpc 2 (x,y) = Y d k e lkx xl pc2 (y) + 

k =1 
-boo 

£ 4e- fa xr 2 (»)> (5) 

/c=-A^/qp c 2 “b 1 

with M qpc2 being the number of conducting transverse modes in QPC2. The transverse modes 
Xk for both QPCs were calculated using the method of Ref. 32 . The scattering amplitudes 
bki dk and the amplitudes of the incoming modes ak of Eq. Q once established are used to 
calculate the transmission probabilities. Throughout this paper we keep Ep = 8meV which 
gives the value of Ap = 2n/kp = 2nh/ a/2 m e sE F ~ 53nm. We choose discretization grid spacing 
Ax = 4nm, which is small compared to Ap- 

After solution of Eq. [2] for each incoming mode the source-drain conductance is evaluated 
from the Landauer formula 

M q pd 

G = Go Y 

i= 1 


( 6 ) 
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where T* is the transmission probability of the i-th mode incoming from the QPC1 to the QPC2 
and Gq = 2 e 2 /h. We refer to the sum of backscattering probabilities as the ’’resistance”, which 
is given by the formula 

-Mqpcl 

R = G 0 J2 R i> ( 7 ) 

2=1 

where Ri is the backscattering probability of i-th incoming mode to the QPC1. 

We discuss below the maps of backscattering probability (SGM-R) and conductance maps 
(SGM-G) for current flow between QPC1 and QPC2. In order to calculate the SGM-R/G 
images we evaluated the ’’resistance” R and conductance G by scanning the system with the tip 
potential given by Eq. ([3]) inside the rectangle shown in Fig. [l](b) . The spacing between two 
subsequent positions of the tip was 4nm, thus in order to obtain the SGM image, presented 
later in paper, for the scan area 400nm x 600nm one has to solve the scattering problem 15000 
times. 


III. RESULTS AND DISCUSSION 

A. In the absence of the tip: Young interference and interference due to QPC2 

detector 

We begin the discussion by setting equal widths of both input and output channels within 
QPCs, u’qpd = '?Uqpc 2 = 40nm. In both the input and the output channels we have a single 
subband transport for the applied Fermi energy of 8 meV. According to our previous paper^ 14 
the single subband transport on the input part allows for resolution of the Young interference 
pattern by the SGM technique. 

The amplitude of the scattering wave function in the absence of the tip is shown in Fig. 
0 b )- Besides the well-resolved Young interference pattern with five more or less straight lines 
of constructive and destructive interference for the waves passing through both the slits [see 
Fig-i c)] one also notices additional vertical interference fringes. This feature is caused by the 
presence of the output QPC2 which reflects the incoming wave function back to the double slit 
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device. 

The role of QPC2 in the scattering is more clearly visible in the real part of the wave 
function displayed as Fig{l|c). We found that a very similar wave function can be reproduced 
by a superposition of three point sources positioned at the center of each slit of the system 
following the Huygens-Fresnel principle for the circular waves propagating from the input slits 
and QPC2 - the latter as due to backscattering of the incoming wave 10 . In the limit of thin 
slit the Fermi level wave function at position r traveling from the slit located at position r source 
can be well described by the angle-modulated Henkel function 


t 'Herikel (h ^source; ^source) 

giAlp |r r S ource | 

^source cos(6*)- 


\/k F \r - 


( 8 ) 


where 6 is angle between r — r source vector and the x axis, i.e. cos(6 ) ) — \x — x SOUTCe \/r and a source 
is a scattering amplitude. For the superposition of the three point sources one has 


~ ^Henkel(T'Wl, 1) + ^Henkel^, r 2 , 1) + 

^Henkel (j’t ^3; ^reflection) j (9) 

where rq = (llOnm, — 80nm), r 2 = (llOnm, +80nm) and r 3 = (670nm, Onm) are the positions 
of slits in the system of Fig. 00- We set a re flection = «source/2, for which we get the best 
agreement with the exact solution obtained from Eq. ([2]). The real part of the wave function 
obtained from Eq.(j9]) is plotted in Fig. [ljd) . Note that the fitted value of a re fl e ction is quite 
large, which suggests that the backscattering from QPC2 will have a substantial influence on 
the SGM images discussed below. 

The Young interference [Fig. [2^d)] and the interference due to backscattering by the QPC2 
detector [Fig. 0b- c)] are present already without the tip. The Young interference involves 
superposition of waves passing through each of the slits with the modulation of the scattering 
density given by 


p = cos(fci7’(r / 1 — r' 2 )). 


(10) 
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B. Interference mechanisms due to the tip 


The backscattering by the tip introduces additional interference effects which will be dis¬ 
cussed below: with the electrons reflected to the source QPC [Fig. [2j[a)] , the tip-induced double 
slit interference [Fig. |2^b)] and the interference of the direct wave with the scattered one [Fig. 
[2j[c)] . The conductance map fringes due to the scattering of the type given by Fig. [2j^a) is well 
known and has been discussed in a number of papers 5 12 . The backscattering (R) is enhanced 
when the phase shift of the wave going to the tip and back [red arrows in Fig. [2j[a)] produces 
a constructive interference with the incident wave [blue arrow in Fig. |2](a)] at the input slit. 
The R signal is then proportional to 

/ = cos(fc F (ri+ ri)), (11) 


where in the argument 2r\ is the phase acquired by the electron wave function on its way 
from a slit of QPC1 to the tip and back from the tip to the same slit. The interference of Fig. 
[2](a) leads to the angle-independent modulation of the R maps which oscillate as a function of 
distance from the slit with the period of half the Fermi wavelength A F /5P^3. 

The tip induces a double-slit interference of the wave passing through one of the slits of 
QPC1 and scattered by the tip with the wave incident from the other QPC1 slit [Fig. |2](b)]. 
Then, an enhanced backscattering can be expected when the interference of the incident and 
returning path is positive, or the R signal will be proportional to 


/ = cos (kp (n + r 2 )). 


( 12 ) 


The result of formula (12) is plotted in Fig. |3](g), with a good agreement with the oscillations 
found in the results of the scattering problem of Fig. [3|e-f). The double slit interference 
according to mechanism of Fig. [2](b) and the resulting SGM effects were never discussed 


before. The equations (11 and [12]) produce the same periodicity of Ap/2 at a large distance 
from QPC1. 

The tip induces also interference of the waves incident from one of the QPCs and backscat- 
tered by the tip - as depicted in Fig. [2](c). The modulation of the scattering density for a fixed 
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tip position is given by 

p = cos (k F (r[ - r 3 )) (13) 


are displayed in the lower panel of Fig. [2](c) for a given tip position. This type of interference 
leads to lateral fringe patterns in SGM-G images to be discussed below. 

C. Resistance maps 

In Fig. |3](a-f) we show the SGM maps of resistance (SGM-R) as functions of the position of 
the AFM tip. The area of the scan is shown by the blue rectangle in Fig. [l](a,b) . We consider a 
single or both QPC1 slits open as illustrated schematically in the insets in the top-right corner. 
For each system we plotted in the second row [Fig. [3jdi-m)] the corresponding probability 
current distribution. In Fig. |3](a,b) we show the results of SGM-R images obtained for a single 
(lower) input slit open. In the absence of QPC2 [Fig. [3](b)] we observe circular fringes in the 
SGM map due to the interference of type given by Fig. [2j)a) between the wave incoming from 
the slit and the wave backscattered by the AFM tip [see the arrow in Fig. 3(d), Fig. |2](a), and 
Eq. [II). 

Figure [3](d) presents the sum of SGM-R maps of Fig. |3](a) and (c) obtained for a single slit 
open. The sum is quite different from the image obtained for the system where both the QPC1 
slits are open (see Fig. [3|e)), which is a signature of the double-slit interference effects, and 
involve both the Young interference [Fig. [2j)d)] and the tip-induced interference [Fig. |2](b)]. The 
Young interference is more clearly resolved in the absence of QPC2 [Fig. |3](f)], although it is 
also present when the QPC2 detector is a part of the setup [Fig. |3](e)]. In order to demonstrate 
this closer we extracted the enlarged fragments of Fig. [3|e,f) in Fig. |4](b,c) . For a simple sum 
of SGM images for separate slits [Fig. [3](d)] instead the elliptic fringes or Young pattern we 
observe a checkerboard pattern [Fig. |3](d) and Fig. 0 a)]. This checkerboard pattern should 
be observed in the experiments for a large number of incident subbands, i.e. at higher Fermi 
energy or wider input channel as discussed in Ref. 11 . 


with r 3 = rq + r tip _ r ^, and r tip _ r ^ is the distance between the tip and r\. The results of Eq. 


13 



11 


a) 



I = cos(/c F (ri 4- ri)) 
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FIG. 2: Schematic presentation of the discussed interference processes. Upper panels show the interfer¬ 
ing paths (blue and red), and lower panels plot the cosine of the phase difference as a function of the tip 
position (a,b) for the SGM-R maps (a,b) or - in the absence of the tip (c,d) - for the scattering density. 
The following interference mechanisms are discussed below: a) due to the scattering of the incoming 
wave on the tip potential wave function may be reflected backward the input slit. Such process will 
lead to circular fringes around the QPC1 slit entrance in the SGM-R images - see the lower panel and 


Eq. 11 b) tip induced double-slit interference when the electron wave is backscattered by the tip to 


the other slit that leads to elliptic fringes visible in SGM-R images - see the lower panel and Eq. (12) 
c) tip-induced interference of unperturbed wave with the wave reflected from the potential tip (this 
process leads to the lateral fringe pattern in SGM-G images and upon additional backscattering by 


the detector to the lateral fringe pattern in SGM-R images) - see Eq. (13), d) Young-type double slit 
interference of two coherent circular waves (characteristic beams visible in SGM-R images) - see Eq. 


(10). 
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In Figs. [3^h-m) we plotted the current density within the system in the absence of the tip. 
Note, for the double slit systems [Fig. [3^1, m)] that the current distribution is very similar with 
or without QPC2. The deformation of the Young interference for the system with QPC2 is due 
to the lateral interference pattern involving the paths marked in Fig. |2](b) which are introduced 
by QPC2. 


The Young interference pattern given by Eq. (10) and calculated for the Fermi wave vector 
is displayed in Fig. |3](n) with a good agreement with the probability current distribution of Fig. 
[3](l,n) and the features of the SGM map of Fig. |3](f). A resemblance to Fig. |3](e) - for QPC2 
present - can also be spotted, although the lines of flat R in Fig. |3](e) are curved. The curvature 
as well as the presence of the lateral fringes can be reduced by placing the QPC2 further from 
QPC1, thus reducing the backscattering by QPC2 but at the expense of the reduced G map 
contrast (see below). 


D. Source-drain conductance maps 


The SGM maps of source-drain conductance (SGM-G) of Fig. [5] exhibit a valley of minimal 
values along the shortest path from the source channel to the drain. The pattern of the image 
obtained for both slits open [Fig. |5](d)] is exactly the same as the one given by the sum of 
images Fig. [5](a) and (b) (see Fig. [5](c)) . This shows that in SGM-G images the double-slit 
interference is absent, which may be quite surprising, given that the SGM-R image clearly 
resolved the interference. 

In Fig. |5](a) a lateral fringe pattern is observed along the classical trajectory. This valley of 
minimal G values coincides with the line of maximal backscattering R observed in Fig. 3(a). 
Figure [6](a) shows that the tip separates the electron wave into beams, which arise from the 
interference of waves incoming from QPC1 and scattered by the tip - as illustrated by Fig. 
[2j)c) . The variation of the electron density can be approximated by Eq. 13 The results of 
formula ( [l3| ) are plotted in Fig. |6](b) with a very good agreement with the numerical result 
of Fig. [6](a), which however contains also the fringes due to backscattering by QPC2. Figure 
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FIG. 3: First row - (a-f) - shows the SGM-R images obtained by scanning the rectangle area indicated 
in Fig. [ljb) by AFM tip described by potential Eq.([3|. The insets in each figure show schematically the 
device for which the image was calculated. Second row - (h-m) - shows the amplitude of probability 
current distribution calculated for system without external potential of the tip. bottom slit open (a,h) 
and (b,i), upper slit open (c,j), a sum (a,c) and (h,j) is given in (d) and (k), respectively. (e,l) and 
(f,m) show the results for both input slits open. In (e,l) QPC2 is present, and in (f,m) it is removed. 


(g,n) are the images obtained from Eq. (10) and Eq. (12). The results were obtained for Uf ip = 15 
meV and dti p = lOnm. 


[6j^a) corresponds to the position of the tip above the point A and explains the valley of low 
conductance in the SGM-G image (e.g. see Fig. |5^a)), since QPC2 is in the shadow cast by 
the tip. Figures |6](c-d) correspond to tip above points B and C for a single slit and explain 
the increase in the conductance around the classical path, which is well visible in Fig. [5|a) 
and (b). If we move the tip to positions B or C (see Fig. [6](b) and (c)) the electron wave is 
still separated into two main beams with one of these beams reaching QPC2, which produces 
the characteristic envelope of high conductance around the “classical path”. The minima of 
the lateral pattern of Fig. [5] appear when a nodal line of the interference passes to the QPC2 


detector. 
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FIG. 4: The zooms of Figures [3^ d-f) denoted by blue rectangles. 
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FIG. 5: The SGM conductance maps obtained for only the bottom (a) or upper (b) slit open; (c) sum 
of (a) and (b); (d) map for both slits open. Points A-E marked by the crosses show the tip positions 
for which we plot the probability current distribution in Figjb] 

The points A, B,C in Fig. |3](a) are the same as in Fig. |5](a) and correspond to plots Fig. 
[6|a-d). Note that point B corresponds to a minimum of R and maximum of G , while point C 
corresponds to a maximum of both quantities. The lateral pattern in the R map appears only 
in presence of QPC2 detector [cf. Fig. [3|a) and Fig. §»], and is a result of scattering first by 
the tip and next by QPC2. The effect of the interference of type Fig. [2^b) depends on the 1) 
position of the tip with respect to QPC2 - QCP1 - return path for the backscattered electrons 
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and 2) the incidence angle of the waves deflected by the tip on the edges of the QPC2 electron. 

Let us now consider both slits open with the AFM tip at position where the electron flow 
from the upper slit is totally blocked (see Fig. |6|e)). The current flux from the lower slit will 
be the only one reaching the QPC2. The tip in this manner turns off one of the slits of the 
source channel. Using the symmetry arguments for the considered device we will get the same 
conclusion with second slit blocked and the first transmitting current. Such argumentation can 
explain why at certain points in the obtained SGM-G maps there is no visible interference. 
Nevertheless, there are points where the potential of the tip does not totally block the current 
from any slit, so an interference could be expected. Such a case is presented in Fig. [6](f ), with the 
tip located at a position where the current distribution is almost the same as in the unperturbed 
system (for comparison see Fig. It). We conclude that the double-slit interference is clearly 
visible in the current distribution but not in the SGM-G images. 
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FIG. 6: The probability current obtained for five different positions of the AFM tip denoted by white 
dots and marked in Fig. [5] (a,c,d) Result for bottom slit open obtained for points A, B and C in Fig. 
[5](a). b) result obtained from simple model given by Eq. (13). e-f) Results for both slits open obtained 
for points D and E in Figj5](d), respectively. 


The lack of Young interference in the source-drain conductance images can be explained 
using a reversed bias and the current flowing in the reverse direction, i.e. from QPC2 to QPC1. 
The upper row of Fig. [7] shows the SGM-G, SGM-R and probability density for the current 
flowing from the double slit to QPC2, while the lower presents the quantities for the opposite 
current direction. We can see that the results for SGM-G are exactly identical for both current 
































16 


directions [Fig. |T](a,d)] in spite of the fact that the electron scattering is very different in both 
setups with a clear Young interference in Fig. [Tj[c) with no similar counterpart in Fig. [T^f ). 

The SGM-G images for both the cases [Fig. [T][a,d)] is bound to be identical due to the 
Onsager microreversibility relation for the single subband transport Tq PC 2,qpci — Tqpc\,qpC2, 
which in terms of the Landauer approach implies no current flow for zero bias. For the electron 
incident from the right there is no reason to expect a presence of the Young interference, which 
is indeed missing in Fig. [T^f ). The absence of the double slit interference in SGM-G follows 
from this observation and the microreversibility relation. Note, that SGM-R image for the 
reversed current orientation [Fig. [T][e)] clearly indicates the classical paths that the electron 
can follow from QPC2 to QPC1. Note, that QPC2 is further 100 nm to the right of the end of 
the figure and that two bright beams are emitted from QPC2 which is only 40 nm wide. 

There is a way to restore the interference in the SGM-G images. Each of the M; n subbands is 
with a certain probability backscattered, transferred to QPC2 or goes out of the computational 
box through the transparent boundary conditions to the rest of the system [see Eq. ®l and 
M in is independent of the tip position. Now if we consider that w qpc2 becomes wider, the ratio 
Grest/G decreases, since QPC2 will be able to transfer more current and the probability of 
the electron transfer form QPC1 to QPC2 increases with the width of the latter. Thus for a 
large width w qpc 2 the value of G res t in Eq. |l| will be small and hence ^ M in ps R + G. We 
can express the conductance of the system in as G ~ ^ M in — R oc — R. We can see that 
for large values of w qpc2 the SGM-G should start to exhibit the interference pattern as the 
SGM-R images do. In Fig. [8]( a-e) we show the results obtained for a large value of QPC2 width 
Wqpc2 = 800nm. Fig. ga) shows the SGM-G image for the bottom slit open, and Fig|8](b) is 
a sum of the images obtained for bottom and upper slit open separately. Now, the SGM-G 
image for both slits open (Fig. |8](c)) is clearly different from Fig. [8](b) with distinct stripes due 
to the Young interference. We can compare this result with SGM-R image of Figjsjhl) which 
is now quite similar to the SGM-G image, particularly in the center and close to QPC2. In 
this region the approximated relation G oc —R is well visible and SGM-G image is indeed a 
negative of SGM-R image. To illustrate this further in Fig. [9] we plotted cross section of these 
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FIG. 7: First row a) The SGM-G image, b) SGM-R image c) and the electron density |^| 2 obtained 
for current flow to the right direction in the absence of the tip. Second row d-f) The same as in the 
first row but with current flow in the opposite direction. The results were obtained for {/tip = 15 meV 
and dti p = lOnrn. 

images along the axis of the device. 

Note that the increased width of QPC2 allows us to restore the interference in the SGM-G 
images, but at the expense of the lost information on the electron trajectories from QPC2 
to QPC1. The double-slit interference is present as long as one does not interfere with the 
measurement trying to determine through which slit the particle passes 1 ®. Here, if we set the 
detector (QCP2) for mapping the electron trajectories - with a small value of ta qpC 2 - we gain 
the information about electron trajectories, but we lose the interference pattern. Increasing the 
width Wqpc 2 leads to reduction of spatial resolution of the detector, so we lose the paths in the 
images but restore the interference pattern. 
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FIG. 8: SGM images for the case of big value of u > qpC 2 = 800nm. a) SGM-G image obtained for bottom 
slit open, b) Sum of images (a) and image obtained for upper slit open (not presented here) c) The 
SGM-G d) and SGM-R image obtained for both slits open, e) electron density \ijj\ 2 obtained for the 
system. Note the similarity of (c), (d) and (e) images. The cross section of (c) and (d) along the blue 
line is plotted in Fig. [9] 



FIG. 9: Cross section of Fig. [8] (c) and (d) taken along the blue line. The black curve shows the 
backscattering R (left axis), and the red one - the conductance G (right axis). 

E. Stability of classical trajectories in SGM-G images 

The classical trajectories as extracted from the SGM-G images should be stable against the 
geometrical parameters of the system, the distance between the QPC1 slits d qpcl in particular. 
For each value of the interslit distance d qpc i we calculated the SGM-G(d qpc i) and SGM-R(d 9pcl ) 
images. The calculations were performed for tu qpc i = w qpC 2 = 40nm. In Fig. [Io|(a) we plotted 
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the Pearson correlation coefficient 22 between SGM-G(d qpc i) and last image SGM-G(400nm) - 
the r(G) curve, and correlation between SGM-R(d qpcl ) and SGM-R(400nm) - the r(R) curve. 
We chose the SGM image for d qpcl = 400nm as a reference one. One may see in Fig. [Toj^a) 
that both lines r(G) and r(R) quickly stabilize around 1 after d qpcl = 150nm, which means 
that through all the values from around d qpc i = 150nm to 400nm both images stay almost 
unchanged. Note, that the double-slit interference in SGM-R images is only well resolved for 
the lowest subband transport in the incident channel 1 -' 1 . 

We found that the SGM-G images are generally stable in function of d qpcl , which means 
that paths are always visible, but there is one exception when the distance between slits is 
small enough such that both trajectories (from each QPC) overlap which makes the calculated 
SGM-G maps difficult to interpret. Sample images of SGM-G and SGM-R for small value of 
d qpc i = 160nm are displayed in Fig. [To|(b) and (e), respectively. For larger values of d qpc i the 
classical trajectories are restored (see Fig. [To|(c) and (d)). Note that for Fig. [lQ|(c) and (d) the 
difference between the values of d qpcl in each case is equal to 16 nm. The results for SGM-G 
are nearly identical. However, the SGM-R images - sensitive to the interference effects - very 
strongly depend on a specific value of d qpc i - cf. Fig. [To|(f ) and (g). The wave function passing 
through both the input slits interferes with the AFM tip, as well as with QPC2. The result of 
the interference in terms of the backscattering depends on the variation of the distance between 
the slits of the order of the period of the waves formed by interference at the Fermi level which 
is equal to Xf/2. 


F. Wide tip potential 

Let us return to the single subband transport in the input and the output leads and consider 
the wider tip potential d tip = 50 nm. The results for conductance G for a single slit - Fig. 
[TT|(a) and both the slits Fig. [TT|(b) still indicate the classical current paths, although naturally 
the width of the G minima is significantly increased. The G map pattern is still very similar to 


the one obtained by a sum of maps for separate slits [Fig. 11 - indicating a lack of double slit 
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FIG. 10: a) The Pearson correlation coefficient calculated between SGM-G/R (r(G) and r(R)) images 
obtained for different values of d qpc i the interslit distance [see Fig. 1, the width of the QPC1 slits is kept 
at 40 nm]. In the inset we show the SGM-G and SGM-R reference images obtained for d qpc i = 400nm. 
b-d) SGM-G images obtained for different distances between the slits, see the label inside each figure, 
e-g) Same as (b-d) but SGM-R images. The results were obtained for t/ t i P = 15 meV and dti p = lOnm. 


interference features in the source-drain conductance maps for M out = 1, as discussed above. 

The ’’resistance” map [Fig. |IT|e)] for both QPC1 slits open contains i) the circular fringes 
near the input slits [single slit interference of Fig. |2](a)], ii) elliptical fringes [double slit inter¬ 
ference of Fig. [2](d)], iii) the lateral fringes [Fig. [2](b)] . All the effects (i-iii) are tip-related and 
in this case dominates the Young interference - which is weaker although still detectable. The 
Young interference is restored when the QPC2 detector is placed further from the QPC1 [see 
Fig. |TIj(h)], but naturally at the expense of the amplitude of the G signal [Fig. [TT|(g)] . The 
reduction of the lateral fringe pattern is observed for R - since it involves backscattering by 
QPC2 - is also observed [Fig. [Tljh-i)] . 
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FIG. 11: The SGM images obtained with U ti p = 8meV and dtip = 50nm. (a,d) SGM-G and SGM-R 
images obtained for lower slit open (b,e) the same but for both slit open and (c,f) images obtained as 
a sum. (g-i) Are the SGM images obtained in situation when the distance between QPC1 and QPC2 
was increased by 400 nm to 1 /iin. Panel (h) corresponds to backscattering in for both QPC1 slits 
open, while in (i) a sum of signals for a single QPC1 slit open. 

IV. CONCLUSIONS 

We have considered imaging of electron trajectories for the double-slit experiment with the 
SGM technique. Several interference mechanisms induced by the tip and the drain contact 
as the electron detector have been found and the paths leading to the interference have been 
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identified. We studied the SGM source-drain conductance maps and demonstrated that the 
classical electron paths are clearly resolved but only for a narrow drain contact, for which the 
double-slit interference features are absent. The double-slit interference pattern is present in 
the conductance maps but only for a wider drain contact, when the electron paths are no longer 
resolved. 

We have indicated that a way to observe both trajectories and interference pattern is to 
look simultaneously at two different SGM maps - for the backscattering R - which contain 
double-slit interference signal and for the conductance G - that reveals the paths. The latter 
allows one to map all the equivalent classical trajectories, but without indicating the specific 
path the electron took on its way to the drain channel. 
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